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ABSTRACT 

The effective-index method and Marcatili’s technique were utilized independently to calculate 
the electric field profile of a rib channel waveguide. Using the electric field profile calculated from 
each method, the theoretical coupling efficiency between a single-mode optical fiber and a rib 
waveguide was calculated using the overlap integral. Perfect alignment was assumed and the 
coupling efficiency calculated. The coupling efficiency calculation was then repeated for a range 
of transverse offsets. 
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1. INTRODUCTION 

The behavior of integrated-optic devices is often polarization sensitive. Polarization-preserving 
optical fibers can be used to couple light with a known polarization into these devices. In this 
study, the theoretical coupling efficiency between a polarization-maintaining elliptical-core fiber 
and a rib waveguide was investigated. To couple light energy efficiently from an optical fiber into 
a channel waveguide, the design of both components should provide for well-matched electric 
field profiles. Due to the complex geometry of a rib waveguide, there are no analytical solutions 
to the wave equation for the guided modes. Approximation techniques must be utilized to de- 
termine the propagation constants and field patterns of the guide. The two methods used here 
are the effective-index method 1-4 and Marcatili’s approximation. 5-7 The refractive index profile 
of the elliptical-core optical fiber was approximated by that of a rectangular channel guide, 8-11 
for which an analytical solution can be determined. 

There are several causes of coupling loss between a fiber and a waveguide. 12-14 These include: 
transverse offset (x,y offset), mode field mismatch, longitudinal separation ( z offset), tilt, and 
reflection. The loss mechanisms addressed in this study are the transverse offset and mode 
field mismatch. Mode field mismatch contributes significantly to interconnect loss, 15,16 hence 
modeling is a necessary tool to design waveguides with low coupling loss. As misalignment in 
the x and y directions is inevitable in practice, the theoretical coupling efficiency results for x 
and y offsets are of interest. 

The integrated-optic rib waveguide and elliptical- core optical fiber that are the subjects of 
this study are shown in Fig. 1 . The single-mode channel waveguide is part of an integrated-optic 



pressure sensor which is based on a strain sensitive ring resonator combined with a micromachined 
silicon diaphragm. This sensor can be interrogated with either quasi-TE or quasi-TM polarized 
light. 17 In this investigation, only the TE mode was analyzed. A single-mode polarization- 
maintaining fiber with an elliptical core 18 was used to excite the sensor at the 830 nm wavelength. 



Figure 1: Rib channel waveguide and elliptical-core fiber. 


2. RIB ELECTRIC FIELD CALCULATION - TWO TECHNIQUES 


For the channel guide, two methods of calculating the propagation constant and electric field 
profiles were utilized: the effective-index method and Marcatili’s technique. The TE modes are 
assumed excited and propagating in the z direction, thus the electric field is of the form 

E{x,y,z) = x ,y) <**-**\ (1) 

where j3 0 is the propagation constant. 


Both techniques provide solutions to the scalar wave equation, 


In Eq. 2, n(x,y) is the index of refraction and k 0 =j^, where X a is the wavelength in free space. 
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2.1 Effective-index method 


The effective-index method reduces the two-dimensional scalar wave equation into two one- 
dimensional problems. First, the rib guide is divided into three separate planar guides oriented 
perpendicular to the y axis, as shown in Fig. 2. Slabs A and C, which have identical charac- 
teristics, correspond to the etched portion of the waveguide, while Slab B corresponds to the 
unetched portion. Ni is the effective index of Slab B and Nn is the effective index of Slabs A and 
C, where N I>It = 0i,n/k o . To determine the degree of confinement in the x direction, a symmetric 
slab waveguide, Slab D, is constructed which is perpendicular to the x axis. Slab D uses Ni and 
N n as the core and cladding indices, respectively, and the rib width W as the core thickness. 
The analytic solution of Slab D yields the effective index, N^ective, of the rib guide. 
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Figure 2: Effective-index method 


The electric field profile of the rib is the product of the field distribution of Slab B, X(x), 
and that of D, Y (y), i.e. 

*(x,y) = X(x)Y(y). (3) 

The effective-index method has intuitive appeal. However, because the guide is analyzed as three 
separate slabs, the field continuity conditions at the interfaces are not met. The results of this 
method are most accurate for ribs with small height and large width. 19 This condition is satisfied 
in this study as the guide etch depth is .09 pm which is small when compared to its 2.7 pm 
width. 


2.2 Marcatili’s method 


For modes far from cutoff, the majority of the field energy is confined in the central core 
region. For purposes of analysis, a pseudo-rectangular waveguide is utilized to approximate the 
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rib structure. 11 The cross section of the waveguide is divided into nine regions as shown in Fig. 
3, in which t is the etch depth and r is the thickness of the etched silicon oxynitride (SiON) layer, 
such that the rib height is r + t. 



Figure 3: Marcatili’s method: pseudo-guide used to approximate rib guide. 


The pseudo-guide is assigned a refractive index profile which approximates that of the actual 
rib. This profile provides readier analysis because n 2 (x, y) is separable in the x and y coordinates; 
it is defined as 11 

n 2 (x,y) = n^ix) + n ,,2 (y ) - n 2 , 


where 


and 


n'(x) = rii for |x| < 


W 

2 , 


n'(x) — 712 


W 

for |x| > — 


( 4 ) 

( 5 ) 

( 6 ) 


n "(y) — n i for — t < y < r, 
n "{y) - for (y < -t) or (y > r). 


( 7 ) 

( 8 ) 


Analogous to the effective-index method, this method reduces the two-dimensional wave 
equation to two separate one-dimensional problems. Marcatili’s technique is considered an ap- 
proximate analytical approach and is most accurate when the waveguide mode is far from cutoff. 5 
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(9) 


Substituting Eqs. 3 and 4 into Eq. 2 produces 


where 


| ^L+kln 2 (x)-fi 

+ k ° n ” 2 M ~ Pi 


X(x) = 0, 

Y{y) = 0, 


Pi + PI ~ k 2 0 n\ = ft. 


( 10 ) 

( 11 ) 


When the mode is well-guided, the majority of the field energy is in the central core (ni region) 
which is a key assumption of this theory. 5 The refractive index profile of the rib waveguide said 
pseudo-guide differ in the outer sha ded region s and sections D and H of Fig. 3. The index of 
refraction in the shaded regions is y2n^ — nf. In regions D and H the actual index is n 2 when 
-t < y < 0 and n\ when 0 < y < r, but in this approximation it is assumed that n 2 is the index 
throughout regions D and H. However, because the modal power in these outer regions is small, 
the mode field profile of the pseudo-guide should closely resemble that of the rib guide. This 
assumption is generally applicable for large etch depths (t >> r). However, for this channel 
guide the etch depth (t) is small, but because the rib is wide (W >> r + t), it is expected that 
the pseudo-guide will provide an accurate field profile. 

Following the development of Varshney and Kumar, 11 the x dependent solutions are 


X(x) = A) cos for M < y, (12) 

X(x) = /4 2 exp for 1*1 > -J" ( 1 3 ) 


Mi = 

M2 = y -*>>?, (15) 


and 6 = 0 for a mode symmetric (antisymmetric) in x. The y dependent solutions are given 

by 


Y(y) = (— jp) 

Y(y) = A c cos (^j+B c sin 
Y(y) = A a exp ^ 


for y > r, 
for — t < y < r, 
for y < -t, 
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( 16 ) 

(17) 

(18) 



where 


76 = tjffi-klnZ, 

(19) 

*5 
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II 

(20) 

£ 
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£ 

(21) 


Assuming quasi-TE polarized incident light, the boundary conditions, 

2 * / , , d$(x,y) W 

n $(x y y) and — continuous at x = ±— , 

ox 2 

$$(x v) 

$(x, y) and — )— — continous at y = —t, r, 
dy ' 


( 22 ) 

(23) 


are applied and eigenvalue equations determined for both the x and y dependent field functions. 
These eigenvalue equations are 


fa' 1 =N-* + (p-l)^ = 0, 

H-2 Ml 


>! 


tan 


-l 


To 

iTc. 


4- tan 


-l 


76 

.7c. 


- 7< (i + “) + (9~ l)* = 0, 


(24) 

(25) 


for the x and y dependent field functions, respectively. For the fundamental mode, p — q — 1 in 
Eqs. 24 and 25. Constants Ai, A 2 , A OJ A&, A c , and B c are determined by satisfying the boundary 
conditions given in Eqs. 22 and 23. 


3. FIBER ELECTRIC FIELD CALCULATION 


The fiber in Fig. 1 was a single-mode polarization-maintaining elliptical-core fiber. Thus, 
the standard Gaussian approximation used for circular-core fibers is not valid for this case. 
Recently, several papers have discussed solutions for elliptical-core fibers. 8-10 In this study, the 
elliptical core was approximated by a rectangle for which an analytical solution can be found. 
The rectangular core dimensions are chosen such that it and the elliptical core have the same 
area and aspect ratio. 9 Therefore, the rectangle dimensions, a and 6 , are given by 

a = Or^- and b' = b (26) 

A diagram of the elliptical and rectangular cores is shown in Fig. 4. Using this rectangle as the 
pseudo-guide, Marcatili’s method (Section 2.2) was used to determine the fiber field profile. The 
field is given by Eqs. 12-21, where t and r have been replaced by b , W has been replaced by 2a , 
and B c = 0. 
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4. COUPLING EFFICIENCY CALCULATION 


Once the electric fields of the rib guide and fiber were calculated, the efficiency of coupling 
between the two components was determined. Because the index of refraction difference between 
the fiber and waveguide was small, Fresnel reflection loss was neglected. 15 The coupling efficiency, 
rj, was calculated using the overlap integral, 14,17 

[jJ$/foy)Sg(g,y) dxd v\ ^ 

1 [/ / $/(*, y)®}(*> y ) dxd y] [/ / y)$z( x > y) dxd v ] ’ 

where the integrals are taken over all space and $f(x,y) and y) are the electric field 

profiles of the fiber and rib waveguide, respectively. The major axis of the elliptical-core fiber 
was oriented parallel to the x axis of the guide. The position of zero x and y offset is chosen to 
correspond to alignment of the maximum field positions of the fiber and rib. 19 After solving for 
the aligned case, various x offsets were introduced and a new coupling efficiency calculated for 
each. This procedure was then repeated for various y offsets while the x offset was zero. 

5. RESULTS 


Field results for the elliptical fiber field using Marcatili’s method are discussed in section 5.1. 
For the rib guide, the field profiles calculated using the effective-index method and Marcatili’s 
approximation are described in section 5.2. The results of the theoretical coupling efficiency 
calculations for various x and y offsets are described in section 5.3. 
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5.1 Waveguide field results 


Theoretical field results using the effective-index method and Marcatili’s approximation are 
shown in Figs. 5 and 6. Results using the effective-index method tend to overestimate the width 
of the field in the x direction 5 as illustrated in Fig. 5. 



Figure 5: Effective-index method: electric field contours for rib guide. 



Figure 6: Marcatili’s method: electric field contours for rib guide. 
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5.2 Fiber field results 


Utilizing Marcatili’s method, the fiber field profile was determined using the parameters n^e 
= 1.484 and n dad = 1.450. The major and minor axes (2a and 2b) of the elliptical-core fiber 
were 2.6 fim and 1.3 /im, respectively. A contour plot of the electric field is shown in Fig. 7. 
For all the contour plots in this paper, the contours indicate regions of equal field magnitude in 
increments of 0.1, with the maximum field magnitude being equal to 1.0. 



Figure 7: Fiber electric field contour plot with elliptical fiber core superimposed (dashed 

ellipse). 


5.3 Coupling efficiency results 


Once the alignment position was determined and the coupling efficiency at this position was 
calculated, an offset in the x direction was introduced. This offset ranged from -3 /x m to +3 /x m 
in .20 /im increments. The coupling efficiency as a function of x offset is shown in Fig. 8. 
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COUPLING EFFICIENCY vs. X OFFSET 



Figure 8: Coupling efficiency vs. x offset. 


COUPLING EFFICIENCY vs. Y OFFSET 



Figure 9: Coupling efficiency vs. y offset. 

Similarly for the y direction, keeping the x alignment, the fiber position was varied from -3 
/i m to +3 /im from the position of alignment. The theoretical coupling efficiency results vs. y 
offset are shown in Fig. 9. 

From the graphs in Figs. 8 and 9, it is apparent that alignment was more critical in the 
y direction. This is expected since the rib field calculated was more tightly confined in the y 
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direction than the lateral direction ( x ) . The maximum theoretical coupling efficiency calculated 
was 81% using the effective-index method and 88% using Marcatili’s technique. The actual 
coupling efficiency is expected to be somewhere between these two results. The difference was 
due to the weaker lateral rib field confinement predicted by the effective-index method, compared 
to the tighter confinement predicted by Marcatili’s method. Because the fiber field profile more 
closely matches the rib field predicted by Marcatili’s method, the theoretical coupling efficiency 
is higher than that calculated using the effective-index method. 

6. CONCLUDING REMARKS 

Theoretical coupling efficiency between an elliptical-core fiber and a rib guide was determined 
for a range of transverse offsets. Two techniques were used to calculate the theoretical electric 
field profile of the rib waveguide: the effective-index method and Marcatili’s technique. For 
this channel structure, both methods produced similar results. The electric field profile of the 
elliptical-core fiber was determined using Marcatili’s method. 

Calculations of this kind are useful in optimizing the fiber choice and rib guide design in 
order to maximize the coupling efficiency between the two components. 
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